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Abstract 

Modeling glueballs as bound states of transverse constituent gluons allows to understand the 
main features of the lattice QCD glueball spectrum. In particular it has been shown in previous 
works that the lightest C-even glueballs can be seen as bound states of two transverse constituent 
gluons interacting via a funnel potential. In the present study we show that such an effective 
potential emerges from the available lattice QCD data. Starting from the scalar glueball mass and 
wave function computed in lattice QCD, we indeed compute the equivalent local potential between 
two transverse constituent gluons in the scalar channel and show that it is compatible with a 
funnel shape, where standard values of the parameters are used and where a negative constant has 
to be added to reproduce the absolute height of the potential. Such a constant could be related to 
instanton-induced effects in glueballs. 
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I. INTRODUCTION 



The study of glueballs is an active domain of hadronic physics, both experimentally 
and theoretically. Much experimental effort has been devoted to the identification of a 
clear glueball signal, but no unambiguous candidate has been found up to now It is 
nevertheless generally accepted that a state with a dominant glueball component should 
be present among the many known f resonances [2|. From a theoretical point of view, a 
remarkable achievement has been the computation of the glueball spectrum thanks to lattice 
QCD calculations I3|, 1^1 . Many works have since been devoted to the understanding of this 

n 

spectrum within different frameworks: Coulomb gauge QCD [5], AdS/QCD duality [6j, 
potential models [3, S| , etc - see also Ref . for a review. 

In potential models, there are two basic ingredients. First, one assumes that glueballs 
are bound states of constituent gluons. Second, one introduces an interaction potential that 
aims to describe at best the QCD interactions. Let us discuss these two points. 

As it is shown in Ref. [h3], the structure of the whole lattice QCD glueball spectrum 
can be understood by assuming that a particular J PC glueball can be labeled by a given 
number of transverse (helicity-1) constituent gluons. The low-lying C = + (— ) states are 
then seen as mainly two- (three-) gluon states. The problem of using longitudinal (spin-1) 
constituent gluons is that light l ±+ states are then present in the spectrum, in disagreement 
with lattice QCD [?], Q]. However, if transverse constituent gluons are used, the spurious 
states disappear j^], and a simple potential model (two-body spinless Salpeter Hamiltonian 
with funnel potential ar — k/t) becomes able to reproduce with accuracy the results of 
lattice QCD in the C = + sector [11]. Even if the transversality is assumed, the value of 
the gluon mass is still a matter of discussion. There are works in the spirit of Refs. 
arguing that a constituent gluon is massless at the level of the confining Hamiltonian (zero 
bare mass) but then gains a constituent mass because of its confinement into a glueball. 
In this picture, the constituent gluon is a posteriori massive and its constituent mass is 
typically riven by ^ \fp^j ~~ notice that the same argument would apply for u and d quarks. 
In Ref. 12] for example, the constituent gluon mass is computed to be around 300 MeV. 
Other studies 8| rather state that a constituent gluon is a priori massive, that is with a 
nonzero bare mass. The underlying idea is roughly that the nonperturbative effects of QCD 
cause a mass term to appear in the gluon propagator. In the present work, motivated by 
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the results of Refs. 
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11], we will only focus on the case where the constituent gluons have 



a vanishing bare mass. 

The potential term remains to be chosen. The potential energy between a heavy quark- 
antiquark pair has been accurately computed in lattice QCD and exhibits a clear funnel 
shape at the dominant order 13j, supplemented with the usual Fermi-Breit relativistic cor- 
rections [14|- The potential energy between two static sources in color octet (static gluons) 
has also been computed in lattice QCD; it is still compatible with a funnel form, for dif- 
ferent values of the parameters because static gluons are color-octet sources rather than 
triplet ones fl3 |. But, constituent gluons should be seen as relativistic particles and one 
can wonder if such a static potential is relevant to describe glueballs. We have recently 
proposed in Ref. 15|] a new method to extract the effective potential between two "phys- 
ical" gluons (in opposition with static ones) from lattice QCD. The starting point is that 
not only glueball masses can be computed on the lattice, but also glueball wave functions, 
especially in the scalar channel 16|, [l7|]. Assuming that the scalar glueball is a bound state 
of two constituent gluons, the inverse problem can be solved, i.e. to compute the effective 
potential corresponding to the lattice QCD glueball mass and wave function from a spinless 
Salpeter-type Hamiltonian. We have given a detailed resolution of this problem in Ref. jl5l. 
but longitudinal constituent gluons were used. Motivated by the results of Refs. 0, 0, 11] 



favoring transverse constituent gluons, we propose in the present work to reconsider the 
calculations of Ref. [ijj] and see whether the effective potential that emerges is compatible 
with a funnel one or not in the case of transverse constituent gluons. 

Our paper is organized as follows. We first recall the lattice QCD data that we have at 
our disposal in Sec. [Til Then we discuss the expected form of a Hamiltonian describing a 
bound state of two transverse gluons in Sec. IHII and we focus on the scalar channel, for which 
reliable lattice data are currently available. In Sec. [IV] we numerically compute the effective 
gluon-gluon potential and comment its compatibility with a standard funnel form and with 
the set of parameters used in Ref. [ill ] . Finally, we draw some conclusions in Sec. [V] 



II. RESULTS FROM LATTICE QCD 

General field-theoretical arguments have been invoked to prove that the lightest glueball 
is a scalar one, i.e. with quantum numbers ++ fis| . This feature has been confirmed by 
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several lattice QCD studies, where the most recent determination of the ++ glueball mass 
is around 1.7 GeV jjj ljj. 

Glueball masses are not the only observables that can be obtained within the framework 
of lattice QCD: Wave functions can also be computed. A first determination of glueball 
radial wave functions in the scalar and tensor channels has been made in Ref. [16] within 
SU(2) lattice QCD simulations, while an SU(3) generalization of these results has been given 
in Ref. 



17| . In this last reference, it is found that 

m 0++ = 1.680 ±0.046 GeV, 



(1) 



in agreement with Ref. 3] . It can be checked in Fig. [T] that the lattice radial glueball wave 
function in the scalar channel is well fitted by the following expression 15j 

R(r) = 2.5 exp —A (r/r ) B 



(2) 



1.472 GeV" 1 , 

0.883 ±0.045, B = 1.028 ±0.132. 



with r 
A 

These results have been confirmed in the more recent work j^. Notice that R(r) has been 
here normalized in such a way that R(0) = 2.5 so that the plot is clearer in Fig. [T], but 
the particular normalization of the wave function is irrelevant for the computations that we 
make in the following. 

As it is argued in Ref. [3], the lattice wave function can be seen as the Bethe-Salpeter 
wave function of a stationary scalar state made of two constituent gluons, in the rest frame 
of the system. It is thus relevant to identify the lattice wave functions of Refs 
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with the Schrodinger-like wave functions of two-gluon bound states. We refer the reader to 
Ref. (Hi for a detailed discussion of that point. 



III. HAMILTON! AN FOR TWO TRANSVERSE GLUONS 

The problem we want to solve in this work is to find a Hamiltonian H such that 

H R(r) = m ++ R(r), (3) 

where R(r) and mo++ come from lattice QCD, as recalled in the previous section. Since 
there are strong evidences supporting the identification of the lowest-lying scalar gluleball 
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as a mainly two-gluon bound state [3, [s[ 
spinless Salpeter form, i.e. 

H = 2 ^/p* + V(r). 



11], we choose for H a standard two-body 



(4) 



We assume the constituent gluons to be transverse and massless as in Ref. 11[ and the 
potential V(r) to be local and radial. Note also that p 2 = p 2 . + L 2 jr 2 . 

Once the radial quantum number n and the matrix element of the square orbital angular 
momentum (L 2 ) are specified (n = here since we deal with the ground state), the only 
unknown quantity in Eq. ([3]) is the local central potential V(r). As we will show in the 
next section, such a particular inverse problem can be numerically solved with accuracy. 
But, before making explicit calculations, one has to wonder about the value of (L 2 ). In 
our previous work [ijj], we considered, as it is the case in most of the potential models, 
that the constituent gluons were longitudinal. In this case, the lowest-lying scalar glueball 
is a 1 1 *S'o) state in spectroscopic notation, thus with (L 2 ) = and (S 2 ) = 0. However, 
this assumption neglects the fact that constituent gluons should actually be transverse. We 
showed in Ref. [lj| that a proper inclusion of the gluon's helicity can be performed within the 
helicity formalism, first presented in Ref. [2]]]. The results are then considerably improved: 
The spectrum that is found is in better agreement with lattice QCD, and spurious states due 
to longitudinal degrees of freedom disappear. Followin g th e helicity formalism, the unique 
scalar glueball state corresponds to the quantum state [11] 



|0 ++ ) 



l^o) 



5 A)} 



(5) 



for which it is readily found that (L 2 ) = (S 2 ) = 2. Notice that the unique pseudoscalar 
state is simply |0 h ) = — \ 3 Pq) with also (L 2 ) = 2. It is worth pointing out that the 
particular form of © is model- independent: Only the requirement that the total spin J 2 
and J z are good quantum numbers and that the gluons are transverse is needed to build 
the s pin -angular part of the various two-gluon helicity states from the general formalism of 
Ref. [21]. 

In summary, once transverse constituent gluons are assumed, the problem that remains 
to be solved is the inverse problem corresponding to Hamiltonian (j3J) in the scalar channel, 
that is 



H 



p< 



+ V(r) 



(6) 
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One could have thought that, since explicitly spin-dependent terms are only present at 
the level of the relativistic corrections, the transversality of the constituent gluons would 
not strongly affect the features of the system. However, the difference between spin-1 and 
helicity-1 gluons dramatically affects the Hamiltonian of the system even at the dominant 
order, because (L 2 ) is equal to 2 instead of 0. 



IV. NUMERICAL RESULTS 



Several numerical methods exist to compute the equivalent local potential from a given 
wave function and its corresponding energy but, to our knowledge, only the Lagrange mesh 
method is able to make accurately such computations with a semirelativistic kinematics of 



the form yp 2 [22]. We refer the reader to Refs. 22|, |23| for further elements about this 
numerical technique, but nevertheless recall the main points for self-completeness. 

In the case of radial equations, a Lagrange mesh is formed of N mesh points Xi which 
are the zeros of the Laguerre polynomial Ln(x) of degree N (here N = 100). The Lagrange 
basis is then given by a set of N regularized orthonormalized Lagrange functions satisfying 
the condition fj(xi) oc 5ij. These functions can be expressed in terms of L^(x). The matrix 
elements T^- of the operator yp^, depending on (L 2 ), can be accurately computed in this 
basis. Moreover, for a local radial potential, one has = V(hxi)Sij in the Lagrange basis, 
with h an appropriate scale parameter. It is readily checked that this nice feature allows 
to rewrite Eq. as V(hxi) = f [T^, m n ++ , R(hxj)] , thus to solve the inverse problem by 



22]. In our case, both m ++ and R(hxi) are 



computing the potential at the mesh points 
known thanks to the lattice results (pQ) and 

Formally, the numerical determination of V(r) in Eq. ([3]) is identical to the one performed 
in Ref. [15]. The crucial difference however is that (L 2 ) = 2 this time and not zero as in 



Ref. [15|]. The potential obtained by using the optimal values of m ++, A, and B is plotted 
in Fig. [TJ It exhibits a confining long-range part, and a quickly decreasing short-range part. 
The errors on these three parameters [see Eqs. (pQ) and (J3D] allow the "true" potential to be 
located between two extremal curves, within the shaded area of Fig. [TJ It can clearly be 
observed that the optimal potential obtained in Ref. 15], that is with longitudinal gluons, is 
different from the new one even if the error bars are included. Again, it should be stressed 
that the use of transverse gluons considerably affects the physical features of the two-gluon 
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FIG. 1: Plot of the numerically computed optimal potential V(r) from Eq. ([6]) (dashed line). It 
is computed from the optimal wave function ([2]) (dashed-dotted line) fitted on te lattice data of 
Ref. [17( (circles). The errors on the glueball mass m ++ and on the wave function parameters A 
and B actually allow every potential which is located in the gray area. These results are compared 
with the funnel potential (|7|) for the values ([8]) of the parameters (solid line). The optimal potential 
obtained in Ref. [l5(] with longitudinal gluons is also plotted for comparison (dotted line). 

system. 



These numerical results are qualitatively compatible with a funnel potential of the form 

V f {r) =Car- — -D. (7) 

In this expression, as is the strong coupling constant, a is the fundamental quark- ant iquark 
flux tube energy density, and C indicates the scaling of the energy density, which is different 
for a gluon-gluon system or a quark-antiquark pair. Several approaches support the Casimir 
scaling hypothesis, i.e. C = 9/4 [24J. The constant D can be used to fit the height of 
potential ([7]) on the numerically computed optimal potential. 

Instead of fitting the parameters of potential (I7j) on our numerical results, we find more 
interesting to check if they are in agreement with the potential of Ref. [11], that leads to a 
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nice agreement with lattice QCD in the C = + sector. This last potential is such that 

C = 9/4, a = 0.185 GeV 2 , (8) 
a s = 0.45, D = 0.450 GeV. 



We argued in Ref. llj that this particular value of D actually comes from instanton-induced 



Uj. A look at Fig. Q] 



forces in glueballs, that roughly contribute to the mass as — D P Sj :0 
shows that the funnel potential (J7J) with the parameters (jHj) is compatible with the numeri- 
cally computed potential. We notice that the ++ mass we find by solving Hamiltonian (j3J) 
with V(r) = Vf(r) and with the parameters (jHJ) is 1.724 GeV, in agreement with the glueball 
mass (JH). 

V. CONCLUSIONS 

In this work, we have computed the effective potential between two transverse, massless, 
constituent gluons from the mass and wave function of the ++ glueball obtained in lattice 
QCD. We followed the method described in Ref. [ijj], but this time the helicity of the 
constituent gluons has been taken into account. The numerically computed potential is 



compatible with a funnel form with the values of the parameters that we used in Ref. [ll( to 
reproduce at best all the low-lying C = + glueballs (not only the scalar one). The present 
results thus draw a bridge between lattice QCD and potential models, and show how physical 
informations about the relevant effective potential to use can be derived from lattice QCD, 
especially in the light hadron sector. 

A negative constant is needed to lower the funnel potential to the value that is numeri- 
cally computed. Such a negative constant has been argued to come from instanton-induced 



interactions in glueballs in Ref. [HI]. It would be very interesting to have at our disposal 
the glueball wave function in the h channel since in this case the instanton-induced mass 
term is expected to be positive. Thus we should have D w —0.45 GeV in the pseudoscalar 
case, the other parameters being unchanged. Moreover, in the other channels (the tensor 
one in particular), one would expect D pa 0. New accurate lattice QCD simulations could 
be very useful to check this point. 

Finally, let us make some comments about the tensor glueball, for which a wave function 



has also been computed in Refs 



• Q[l2|. 



If constituent gluons were spin-1 particles, the 
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lowest-lying 2 ++ state would be a l 5 ^) one, and its wave function should be qualitatively 
similar to the one of the ++ glueball, since both states would have (L 2 ) = 0. With helicity- 



1 gluons however, (L 2 ) =4 in the tensor channel Hj, and the wave function should thus 
be qualitatively different from the one of the scalar glueball. The data of Ref. are 
unfortunately not accurate enough to make a definitive comparison between the scalar and 
tensor glueball wave functions, although their behavior seems to be rather different. We 
hope that new data will be available in the future, in order to further check the relevance 
of the transverse gluon picture. 
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